NASA Contractor Report 4212 


Numerical Simulation and Comparison 
of Symmetrical/Supercritical 
Airfoils for the Near Tip Region 
of a Helicopter in Forward Flight 


F. F. Badavi 

Planning Research Corporation 
Aerospace Technologies Division 
Hampton, Virginia 


Prepared for 

Langley Research Center 

under Contract NAS 1-18000 


NASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Division 


1989 



ABSTRACT 


A Navier-Stokes (N-S) equations solver named HEFTA (HEIicopter Flow Transonic Analysis) is 
developed to demonstrate the advantage of computational fluid dynamics (CFD) for the prediction of 
airloads and acoustics properties. HEFTA's primary emphasis is on a onebladed helicopter rotor in hover 
or forward flight at transonic tip condition. The solver is based on a three dimensional, time accurate, 
compressible, thin-layer Reynolds-averaged formulation. The equations are solved in an inertial 
coordinate system on a body-conformed curvilinear grid of C-H topology. The major emphasis of the initial 
phase of this work is the study of the aerodynamics of the shock wave boundary layer interaction. Detailed 
boundary layer and global numerical comparisons of NACA-0012 symmetrical and CAST7-158 
supercritical airfoils are made under identical hover and forward flight conditions. The rotor wake effects 
are modeled by applying a correction to the geometrical angle of attack of the rotor blade. This correction 
is obtained by computing the local induced downwash velocity as computed by CAMRAD. The coupling 
of HEFTA and CAMRAD is accomplished in order to place the airfoil in its proper wake environment. The 
calculations are performed on the Numerical Aerodynamics Simulation (NAS) Cray 2 and on the VPS 32 (a 
derivative of Cyber 205 at Langley Research Center) for a model helicopter in hover and forward flight. 
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SYMBOLS 


3oo 

c d’ c dp- c dj 


Cf - 

c p 


speed of sound 

total, pressure and skin friction drag coefficients 

/I 2 

skin friction coefficient, c f =T w/2 p e u e 
pressure coefficient, cp=(p-P«}/ ^-PooU^ 


c 

C 


e 

F ,Fyi G,Gy, H,Hy 


characteristic length scale (chord of the rotor blade, taken as unity) 
constant in Sutherland's equation 

~>2 

total energy per unit volume, nondimensionalized by po ° 800 
inviscid and viscous flux vectors of mass, momentum, and energy 


I 

H 

J 

Moo, M| 

An, AS 
Pr 

P. Poo 


Q 

R 


r 

S 

t 

T 

Uoo. U e 

U, V, w 


identity matrix 

boundary layer shape factor, H = 879* 

Jacobian of coordinate transformation 
free stream and blade tip Mach number 

mesh spacings normal and tangential to the airfoil in physical space 
Prandtl number 

local and free stream static pressures (p nondimensionalized by p °° 800 ) 

conservative variable vector 

radius of rotor blade 

tip Reynolds number R e = p^ujc/poo 

radial distance from center of rotation 

viscous stress vector 

time (nondimensionalized by R/a») 

temperature, also transpose of a matrix quantity 
free stream and edge of viscous layer velocities 
components of velocity in physical space (nondimensionalized by IL) 


U, V, W contravariant components of velocity in computational space 
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inertial and blade fixed coordinates 


X, X b 


x, y,z 

coordinates in physical space 

Y 

specific heat ratio 

aj,« e 

induced and effective angles of attack 

^.c 

transformed coordinates in computational space 

K 

thermal conductivity 

X 

coefficient of bulk viscosity 

p 

viscosity, also advance ratio p = u^QR 

Pe-Pt 

effective and turbulent eddy viscosities 

£x> £y- £z 

metrics of transformation 

P. Pe- P°o 

local, edge and free stream fluid densities 

V 

azimuth angle 

n 

angular velocity of rotor blade 

8* 

boundary layer displacement thickness, 8 = j 1 dz 

O' Pe^e/ 

e* 

boundary layer momentum thickness, 0 = J -^-(1 dz 

J 0 ppj e \ uj 

e c 

geometric angle of attack 

x w 

wall shear stress, x w = 4 hr 1 q 

X 

time in computational space 

x xx> x xy> x yy<- 

viscous derivative terms 

V 

kinematic viscosity 

Subscripts and 

Superscripts 

~,e,w 

conditions at infinity, edge of viscous layer and wall 

i t|, C; x, y, z 

derivatives with respect to n, C: x, y, z 

A 

quantities in computational space 

- 

quantities in rotating coordinate frame 

— > 

denotes dimensional quantities 
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I. INTRODUCTION 


The flow around rotor blades is a complex phenomenon involving transonic flow with associated 
shock wave boundary layer interaction. This results in shock induced flow separation, blade vortex 
interaction on the advancing blade, and dynamic stall on the retreating blade. Other aerodynamic 
complexities, such as rotor fuselage interactions, aeroacoustic effects, and the overall unsteady nature of 
the rotorcraft flow field caused by aeroelastic effects, impose major limitations on the high speed transonic 
flight of rotorcraft. These limitations are manifested through undesirably high vibration levels, power 
divergence, noise, and component fatigue. These different phenomena can not be computed accurately 
by one single method. For this reason, highly specialized methods are required for any particular aspect 
of rotor blade aerodynamics. For many years, helicopter designers have used a mixture of simplified 
integral methods, such as panel and lifting surface computer codes, wind tunnel data, and design charts, 
to perform numerical analyses with the knowledge that the technique of combining traditional integral 
methods, with experimental data for computing rotorcraft aerodynamics loads is unable to compute the full 
impact of transonic effects. 

Numerous investigators have addressed the limitations of integral methods by using nonlinear 
finite difference methods to compute the aerodynamics of rotor blades. Finite difference methods 
provide the technique to compute the entire nonlinear flow field about the rotor blade at transonic 
speeds. The primary difference between linear integral methods and nonlinear finite difference methods 
is that linear solutions depend only on shear layer conditions and blade surface while finite difference 
solutions depend on the entire flow field. The finite difference methods also have the advantage of being 
designed to compute the transonic flow nonlinearities associated with high speed advancing blades. 

They have progressed from the solution of transonic small disturbance equations to full potential 
equations and more recently to the Euler and Navier-Stokes (N-S) equations. 

In this work, an N-S solver for calculating the viscous transonic flow on rotor blades in hover and 
forward flight is described. The time accurate, unsteady, thin layer N-S equations are solved in an inertial 
frame of reference, rather than in a blade fixed rotating reference. This provides a unified scheme which 
has the capacity of handling both hover and forward flight in an identical manner. This approach results in 
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an increase of CPU time, but allows the same overall formulation to be used for both hover and forward 
flight. 

For the calculation of viscous transonic flow about hovering or advancing rotor blades, the N-S 
equations describe shock location, pressure and skin friction drag and vortical motion accurately. 
Calculation of vortical motion for the rotor blade flow field is significant since the field is characterized by 
blade vortex interaction (BVI) and ,in particular, by complex vortical wakes. These must be accounted for 
in order to describe the forces acting on hovering or advancing rotor blades. Figure 1 shows the profile of 
the vortical sheets in hover mode as they pass through the rotor blade. For vortical wakes, the contraction 
of the helical wake as it passes through and below into the rotor disk and the interaction of the tip vortex of 
one rotor blade with the following are significant. The vortex distention and contraction depends on the 
downwash velocity which significantly influences the accuracy of loads prediction on the blade. 

In this work, the rotor wake effects in forward flight are modeled in the form of a correction to the 
geometric angle of attack of the rotor blade.This correction, depending on the flow conditions, varies 
along the blade span and azimuth direction, and reduces the geometric angle of attack e c . So far, most 

hover and forward flight calculations have depended on this vortex induced surface inflow boundary 
condition method which has proven to be farily successful for moderate to high advance ratios where total 
flow is much larger than induced flow. The modeling of wake effects is done by the use of the free wake 
analysis program CAM RAD (ref. 1). This choice of calculation procedure is due to the fact that lifting rotors 
in forward flight have wake systems that are too complex to be modeled in the current N-S equations 
solvers. In addition, there is not a straightforward way to couple the blade trim solution directly into any 
finite difference scheme. With this approach and others (ref. 2), estimates of the downwash velocity Vj 
and corresponding induced angle of attack aj are obtained from free wake analysis calculations. Since aj 

has a negative value, the influence of the vortical wake is to reduce the geometric angle of attack so that 
an effective angle of attack a e = 8 C - otj is obtained. 

In CAM RAD aerodynamic loads are calculated by computing induced angles of attack along the 
rotor span which are the result of the trimmed rotor and wake system. The blade lift coefficients are then 
decided through table lookups of two dimensional airfoil data. The coupling between HEFTA and 
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CAMRAD for a set of specified induced angles of attack is accomplished by replacing the table lookups 
with lift coefficients that are calculated from HEFTA. This process starts by calculating a trimmed solution 
from CAMRAD where the spanwise lift coefficients are completely obtained from airfoil tables. CAMRAD 
then determines a spanwise set of wake induced angles along the rotor surface. The calculated inflow 
angles are then used to model the wake effects in HEFTA. HEFTA then produces a set of unsteady lift 
coefficients which are applied to the subsequent CAMRAD trim solution stage. This iterative process 
continues until the values of induced angles of attack does not change. 

In this work, a previously developed .three-dimensional, fixed wing, thin layer N-S equations 
solver (refs. 3) is modified to simulate the transonic flow field about rotor blade. An inertial coordinate 
system is used on a body conformed curvilinear grids around the blade. A modular solution approach is 
adopted where an independently generated, body fitted grid of C-H topology (C, in wraparound radial 
direction and H, in the spanwise direction) is coupled with the N-S equations solver. The three 
dimensional C-H grid was generated by assembling two dimensional spanwise sectional C-grids. The two 
dimensional C-grid around the blade is generated by the application of an in-house algebraic grid 
generation procedure. The N-S equations solver is developed with an implicit, upwind-biased, finite 
volume scheme which uses the Van Leer's flux vector splitting and the Beam-Warming’s approximate 
factorization of the conservative form of the Reynolds-averaged, thin layer N-S equations to obtain the 
flowfield quantities. Modifications were necessary to implement the rotation matrix into the flow solver to 
accurately solve the rotor blade flow field. 

Part II describes the governing equations and numerical scheme. Part III presents the 
computational results. Part IV describes the near future plans for the present algorithm and the 
conclusions. 
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II. GOVERNING EQUATIONS AND NUMERICAL SCHEME 


Problem Formulation 


The time dependent compressible N-S equations express the conservation of mass, momentum, 
and energy for an ideal gas in the absence of external forces. The strong conservative law form of the N-S 
equations are used for shock capturing purposes. The equations in Cartesian coordinates in 
nondimensional form can be written as 


ao + 3F + dG t 3H _ d p v t dGy f dH v 
at ax ay az ax ay az 

where 

Q = (p, pu, pv, pw, e} T 

/ 2 \ T 

F = |pu,pu +p, puv, puw,u(e + p) f 

I 2 \ T 

G=^pv, puv, pv +p,pvw, v(e + p)^ 

/ 2 l T 

H = | pw, puw, pvw, pw +p,w(e + pw 
F v =Re {o,x X x> T yx- x zx>Px} 

G V = Re {0,T X y,tyy,T z y, Pyj 
Hy*Re {0,t XZ ,Ty Z ,'t zz ,p z | 


( 1 ) 


( 2 ) 


with 
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T xx~ ^(u x + Vy+ w z ) + 2pu x 

X yy = ^U)(+Vy+W 2 j + 2pV y 

"t Z^ = ^|u X + Vy+W Z |+2^W Z 

T x y=Tyx = M|Uy+V x } 

T xz =t zx = M{ u z +w x) (3) 

Xyz=X Z y=M|Vz + w y) 

P x =r KPr 3 x e,+ UT xx +VT xy +WT xz 
-1 

Py=YKPr 3ye|+UXy X + VXyy+WXy Z 

Pz=Y KPr'^zei+UTzx+vxzy+wxzz 

■1 J 2 2 2\ 

e|=ep -0.5\u +v +w / 


The Cartesian velocity components u, v, and w are nondimensionalized by the free stream speed of 
sound a^. 

A general three dimensional curvilinear transformation between the Cartesian variables (x,y,z) and 
the generalized coordinates (^,ti,Q is applied to (1), where C corresponds to the coordinate normal to the 
body surface. The vector Q represents density, momentum, and total energy per unit volume, and p is 
the pressure defined from the equation of state of an ideal gas as 


P 


= (Y-1) 


e - 



2 2 
+ v 



(4) 


The transformation from (x,y,z) to (!;,t|,Q is defined as 


x = t 


% = 4(x,y.z,t) 

T] = Ti(x,y,z,t) 

C = C(x,y.z,t) 
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( 5 ) 



The resulting transformed equations in nondimensional form are written as 


where 

6 = J ^p.pu, pv, pw,e} T 

~ -1 ( \ T 

F = J |pU, puU + £ x p, pvU + S y p, pwU + £ z p, (e + p)U -4 t P} 

G = J ^pV.puV + iixP.pvV + Tiyp, pwV + Ti z p,(e + p)V-Ti t pj (7) 

H = J ^pW.puW + CxP.pvW + CyP.pwW + CzP^e + pJW-Ct pj 


and 

u = ^t+^x u+ V + ^2 w 
V = Tlt + T lx u + T 1y v + T l 2 w 
w =Ct + Cx u + C y v + Cz w 


where U, V, and W are contravariant velocity components written without metric normalization. The 
corresponding viscous flux terms are 

~ -1 -if 

F v = J Re | 0,^ x t xx + ^yT X y + ^ z X xz , ^ x Ty X +^yTyy+^ z Xy Z , ^ x T zx +^yX Z y+^ z t zz , ^ x p x +4yPy+^zPz| 

~ -1 -1/ 

Gy = J Re / 0, lf X X XX +TlyX X y+TJ z X xz , 11 x Xy X + TlyXyy + Tt z Xy Z) T| x x zx + T|yX Z y+ T| ^ zZ > T1 x Px + T1 y Py + T1 Z Pz 

H v = j' 1 Re* 1 |o,C x x xx + Cy T xy + Cz x xz> Cx T yx + Cy x yy + Cz T yz> Cx T zx + Cy t zy + tz T zz> CxPx + CyPy + CzPz j 

( 9 ) 
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where the components of heat flux vectors and shear stress tensors in nondimensional form were given in 

(3). 

Applying chain rule relations, the Cartesian derivatives in generalized coordinates (!j/n.C) become 
u x =^ x U4+Ti x u T1 + CxUr 

Uy=SyU£+11yU T1 +CyU{- (10) 

u z =^ z u^+ti z u t1 + C z uc 

and the transformation metrics become 

^x= J(y^c-ypri) 

5y- J N‘C-*n z c) 

Cx-J^iT^n) 

?y J ( x i z r x ?n) 

C 2 -j(x^-y ?n ) (11) 

ix«J(^c-y?c) 

, 'y= J (x?C-x?5) 

nz-Jfy^-x^ 

^t = _X T^X _ yT^y~ Z T^Z 
^{“-Xx^x-yx^y-ZT^z 
Ct = - x xCx-yx Cy _z xCz 

where J is the Jacobian of the transformation and is given by 
J^x^xC z~ T lzCy) + 4y( T 1zCx -T lxCz) + £z( 1 1xCy _ 'nyCx) 
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Since the dominant viscous etfects of high Reynolds number turbulent flow arises from viscous 


diffusion normal to the body surface, a thin layer assumption is employed. Only the viscous diffusion 
terms normal to the body surface are retained and all viscous derivatives along the body in r\ direction 
are. dropped. That leads to thin layer N-S equations as 

9Q 9F 3G 3H R -13S (13) 

jr + * + 5f * at 

where 

S = j' 1 |o, M (Cx + Cy+Cz)^C +(M/3)(CxUC + CyV C+ C z w C )C x , M^x + Cy + C z ) v; 

+ (m/3)(C x U£+C y V£+C Z W£) Cy, 4^X + Cy+ C z )w^+(M/3) (C x U^+CyV^+C z W^, j(cJ+Cy+C z ) 

x 0.5p(u 2 + v 2 + w 2 )c + pPr' 1 (Y- l)" 1 (afjc + (p/3) (t x u+Cy v + Cz w ) x (CxU£+CyV(;+Cz w c)j 

(14) 

In the above equations, u, v, w are the velocity components in the inertial rotor reference frame 

(x,y,z,t) and (13) is solved in this frame. The inertial coordinate X= (x,y,z,t) is related to the blade fixed 
coordinate Xfo = ( x .y. z .t)by 

X(x,y,z)=R(t)Xb(x,y,z) (15) 

t = t 

where R(t) is the rotation matrix (ref. 4) and is given by 

cos fit -sin fit 0 

sin fit cos fit 0 

0 0 1 
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and Qt represents the azimuthal sweep of the rotor blade. This establishes a relation between metric 
quantities in the inertial and blade fixed frames as (ref. 5) 

^ = £ x cos Qt sin Qt 
£y = £x sin Qt ' + ^yCOS Qt 

Ti x = T^cosQt-TiySinQt 
T)y= iixSin Qt + TiyCos Qt 

T\z =r \z (17) 

C x = Cx cos ^t- Cy sin Qt 
Cy=Cx sinft * + CyCOS Qt 
Cz=Cz 

4t=Qy^x-«x^y 
T| t = Qynx-Qxtty 

Ct=QyCx-fixCy 

Equation (13) is closed by Stokes' hypothesis for bulk viscosity as 

Jl + ip = 0 <1«> 

3 


and Sutherland's law for molecular viscosity as 


JT \ 312 Toc+C 

* \T«) T + C 


(19) 
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where is freestream temperature (T^ = 460 R) and C is the Sutherland’s constant. All calculations are 

done under the assumption that the flow over the entire upper surface of the rotor blade is turbulent, 
which is equivalent to the experimental placement of a transition strip at the leading edge of the blade. 

Spatial Differencing 

/\ XS. /N 

The generalized fluxes F, Gand H , representing pressure and convective terms, are spatially 
differenced by a second order accurate upwind biased flux vector splitting method, i.e., split into 
differenced forward and backward contributions and differenced accordingly. For instance, in one 
dimension the flux difference in the $ direction is (ref. 6) 

8^F=8^F + 8^F [ ' 


where 8^ and 8^ denote general backward and forward divided difference operators. The flux vector in 
three dimensional generalized coordinates is split according to the formulation of Van Leer (refs. 7-8). 

/s 

The flux F_ for example, is split according to the contravariant Mach number in the £ direction as 



with 


( 21 ) 



where V represents the gradient operation. For supersonic flow, where |M^|>1, F is 


13 


^ 

F = F 


F =0 


/N /\ 

F = F 


✓V+ 

F =0 


Me > +1 


Me<— 1 


( 22 ) 


and for subsonic flow where |M^|< 1 , F is 

= ^ {fmasff f TOSs[k4-u±2a4/Y + u] ) 4 ass [k y (-u±2a4/7 + v] , 

^mass[^ z ^ _ u - 2a4/y + w ] . ^energy } 
with 


(23) 


f mass =± P a 4 M ^ ±1 ) 2/4 


± ± 

^energy = ^mass 


(Y - 1 ) u ± 2 (y - 1 ) ua«,+ 2a„ 2 ]/(Y 2 - 1 ) + y(u 2 + v 2 + w 2 ) 


(24) 


The equations above are solved with a finite volume approach. In the E, direction the surface area of the 

cell surface is the cell volume is J' 1 , and (k x . ky. kz) are the direction cosines of the cell interface in the 
J 


£ direction and are defined as 


(k^y-HeS*! 


(25) 
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For = At] = A£ = 1 , the implementation of split flux differences are done by a flux balance across 


the i-th cell holding spatial indices j and k constant as 


— -f- 

5^Fj +8^Fj 





^ (1 +\ 

= 

[f (q )+f (q ) 

i + 1 / 2 - 

[f (q j+F Iq JJ 


i- 1/2 


(26) 


where F (q )i+l/2 refers to a forward flux evaluated using the metric terms at the cell interface i + 1/2, 
and the state variable, Q, is obtained by the upwind second order interpolation of cell centered variables 
as 


Qj + 1/2=3/2Qj-1/2Qj_i 
^i+1/2* 3/2 ^i + 1 ” 1/2 Qi+2 


(27) 


For the diffusion terms representing shear stress and heat transfer effects, the fluxes are 
differenced through a second order central difference scheme where second derivatives are treated as 
differences across cell interfaces of first derivative terms, holding spatial indices i and j constant, as 

/v /\ 

8 C Hy k =Hv k+ 1/2" Hv k- 1/2 (28) 


where 


^ ^ II i 

= ^ 0, <t>iU;+Cx<t>2> 4*1 V (; + Cy < l , 2* < l , l w (; + Cz < l > 2> $1 

l"l 6 [ J 


f 2"L +(Y _1) ^ ^ 


+ W<1>2 , 


(29) 


where q is the heat flux term and 
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( 30 ) 


2 2 2 
<f>1 = Cx+ Cy+ Cz 

<t>2=(CxU;+CyV ; +C z W;)/3 


5 C u k+1/2 =u k+1" u k 


Time Differencing 


The linearized, backward time approximation in delta form for three dimensional equations is given 


by 


I 

JAt 


+ 



+ 




3H V 
3Q , 


AQ = -l(q") 


(31) 


with L(Q n ) as the discrete representation of the spatial derivative terms in (13) evaluated at time step n. 
The above equation (31) is spatially factored (ref. 9) and solved as a series of sweeps through the mesh as 


I . « 3G 
JAt 6n 3Q 


AQ 




I + g 5H __ 3H V 
JaT ? \3Q iQi 


aq**=(— !— Iaq* 

JAt 


-U«J£ 

JAt *dQ 


AQ = — !— AQ 
JAt/ 


(32) 


where 


AQ = Q 


n+ 1 


-Q 


n 


(33) 
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Turbulence Modeling 


Due to the high Reynolds number of the flow, it is assumed that the entire upper surface of the 
rotor blade is turbulent. In the governing equations, the effect of turbulence is accounted for through the 
concept of eddy viscosity. In the momentum equations, the molecular viscosity p is replaced by an 
effective viscosity p e as 


p e =p + p t ^ 

where pj is the turbulent viscosity and the Baldwin Lomax (ref. 10) algebraic turbulence model is used to 
evaluate the turbulent quantities. The turbulent viscosity term, pj, plays an important role in shaping the 
boundary layer and tip vortex formation profiles. In high Reynolds number transonic flow .the rapid 
thickening of the boundary layer influences the location and strength of any shock that may be present. 

Far Field Boundary Conditions 

Depending on the nature of the flow at the boundary (inflow vs. outflow, subsonic vs. 
supersonic) .flow field quantities are either fixed or extrapolated from the interior. A compressible two- 
dimensional line vortex solution (ref. 11) is used for far field velocities that are fixed. This method is only 
used for steady problems due to the time lag of the circulation on the blade and its effect on the far field 
flow. For unsteady cases the far field velocities are set equal to zero when not extrapolated. 

Blade Surface Boundary Conditions 

At the surface of the blade, the no slip conditions are used with the transpiration velocity 
modifications to include the effects of the wake. With the transpiration method the velocity at the blade 
surface is set equal to downwash velocity instead of zero, where downwash velocity is given by 

Vwake = (£ty + V„ sin \|/ cos a TPP ) tan a* (35) 
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Grid Methodology 

A modular solution approach is adopted where an independently generated, body-conformed 
grid of C-H topology (C in wraparound-radial direction and H in the spanwise direction) of 351 x 80 x 21 
cells in the t), £ directions is coupled to the flow solver. This block, with over half a million grid cells, was 
generated by an assembly of 21, two dimensional spanwise sectional C-grids. The two dimensional C-grid 
around the blade is generated by the application of an in-house algebraic grid generation procedure. 
Figure 2 presents the section profiles of the supercritical airfoil CAST7-158 and the symmetrical airfoil 
NACA-0012. The half thickness profiles and camber profiles of each airfoil are also included. Note that 
| the camber profile of the supercritical airfoil is quite different than the straight line profile of the symmetrical 

i 

airfoil. Figure 3 presents the partial view of the C-grid profile for the two airfoils with 351 wraparound and 
80 radial points. In the radial direction the grid lines are heavily stretched to resolve the viscous sublayer. 
The minimum grid spacing in the radial direction is 2 x 10' 5 c (c=chord) which corresponds to the value of 
y + = 3 at a tip Reynolds number of 5 x 10 6 . In the wraparound direction, the important requirements are 

i 

that the grid lines are clustered in regions where steep gradients of flow variables are expected (i.e., at the 

i 

leading edge, trailing edge and near blade tip) and that a smooth and continuous change of the size of the 
adjacent grid cells is maintained. The C-grid is stretched 5 chords from the rotor surface in the wraparound 
and normal direction, and spanwise grid distributions have inner boundaries located 4 chords from the 
rotor tip. The outer spanwise grid boundary is located approximately 2 chords beyond the tip of the rotor 
blade. The rotor blade itself is untapered and untwisted and is not conically faired at the tip (i.e., sharp 
edged tip). In the spanwise direction, a nonuniform distribution of C-grid planes is used to concentrate 
the two dimensional planes near the tip region. A simplified schematic of the entire C-H grid block with the 
embedded blade rotor is shown in figure 4 with C-grid planes extending outward in the spanwise 
direction. 

The C-grid plane is transformed into a generalized coordinate system as 
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y c . z. t) 
C - C (x, y c . z, t) 


( 36 ) 


X = t 


where y c refers to a constant value of spanwise location. The transformation allows the boundary surfaces 
in the physical plane to be mapped into rectangular surfaces in the computational plane. This 
transformation simplifies the boundary condition requirements in the computational domain and allows 
grid points to be clustered in the physical plane in regions of steep gradients. Figure 5 taken from 
reference 12 shows this transformation with the section in the trailing edge enlarged to show the high 
aspect ratio As/An requirements in the viscous sublayer. 

III. COMPUTATIONAL RESULTS 

The primary purpose of this work is to numerically compare the behavior of transonic flow of a 12% 
thick supercritical airfoil section (CAST7-158) with that of a symmetrical airfoil section (NACA-0012) under 
identical forward flight conditions to draw conclusions concerning the blade performance and noise 
propagation in the near tip region of a rotor blade. This is accomplished by performing computations with 
the following flight conditions: 

M t = 0.86 
6 C = 8° 
y = 90° 

H = 0.3 
Re t = 5x 10 6 

y/R = 0.96, 0.90, 0.80, 0.70 


where M t is the tip mach number, 0 C is the collective pitch, y is the azimuthal angle of rotor disk, p is the 
advance ratio, Ref is the Reynolds number at rotor tip, and y/R is the spanwise locations where 

computation is performed. The untapered and untwisted rotor blade has an aspect ratio of 6 and has a 
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CAST7-158 or NACA-0012 airfoil section profile. A typical computation requires 18 x 10 6 words of 
memory and 2 x 10‘ 5 seconds of CPU time per grid point per iteration. 

Performance 

Figure 6 presents the surface pressure distribution for CAST7-158 (solid line) and NACA-0012 
(dashed line) at y = 90° and y/R = 0.96, 0.90, 0.80, 0.70. It shows that CAST7-158 has retained its 
attached turbulent flow profile up to x/c = 0.7 while NACA-0012 has a massively separated turbulent 
profile with a shock position at x/c * 0.4. This difference in shock position is observed at all four locations 
and accounts for the performance of each blade. Table 1 presents the results of the section lift and drag 
coefficients with drag being composed of pressure drag and skin friction drag. Due to a delayed .shock- 
induced, separated boundary layer, the lift coefficient of CAST7-158 is higher than the NACA-0012 at all 
four y/R locations. This accounts for the superior performance of the supercritical airfoil over the 
conventional airfoil at transonic speed. 


Boundary Layer Quantities 

Figures 7(a) - 7(h) present the boundary layer velocity profiles for the CAST7-158 (solid line) and 
the NACA-0012 (dashed line), normalized by the velocity at the edge of the boundary layer, at x/c = 0.01 , 
0.10, 0.44, 0.57, 0.67, 0.75, 0.86, 0.99 and y/R = 0.96, 0.90, 0.80, 0.70. Only the very narrow viscous 
sublayer with a maximum height of z/c = 0.004 is shown to study the behavior of two airfoils very near the 
blade surface. Figures 7(a) - 7(b) show that at x/c = 0.01 and x/c = 0.10 the velocity profiles of the two 
airfoils are somewhat similar. This is due to the fact that near the leading edge of the rotor blade, the flow is 
attached and a separation zone does not exist, which was also observed on the c p plots of Figure 6. 

Figures 7(c) - 7(d) at x/c = 0.44 and x/c = 0.57 are immediately to the aft of the shock position on the 
NACA-0012 and show an unfavorable velocity profile for the NACA-0012 compared to the CAST7-158 for 
all four y/R locations. This unfavorable profile is caused by an adverse pressure gradient behind the shock 
location which causes the u component of velocity to slow down and eventually generate a very weak 
circulation near the surface of the blade as is shown by Figure 7(d) at y/R = 0.96 near the tip. The profiles 
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at different y/R stations at these two x/c locations are also of interest. In figure 7(c) the profiles for the two 
airfoils are somewhat similar near the tip and deviate toward inner locations, while figure 7(d) indicates a 
small region of flow reversal near the tip and a gradual improvement of the NACA-001 2 profile toward the 
inner locations. These differences are due to the fact that the flow over the NACA-001 2 is gradually 
reattaching itself as it moves toward the aft region of the blade. 

A similar situation for the CAST7-158 is shown in figures 7(e) - 7(f) at x/c = 0.67 and x/c = 0.75 with 
a shock position at x/c = 0.7. Here the region of unfavorable velocity profile near the tip is small and 
gradually moves toward the inner locations as shown by figure 7(e), which results in a large reversed flow 
region near the tip at x/c = 0.75 and gradual profile improvements toward the inner locations as shown by 
figure 7(f). Figures 7(g) - 7(h) indicate the possibility of a second smaller reversed flow region for the 
CAST7-158 near the tip as is shown in figure 7(g) with improvements of the flow profiles in figure 7(h). It 
must be noted that at the last four x/c locations, namely x/c = 0.67, 0.75, 0.86, 0.99, the NACA-001 2 has a 
more favorable profile than the CAST7-158 due to the partial reattachment of the flow at these locations. 
Therefore, for the first 60% of the upper surface, the CAST7-158 has a more favorable profile, while for 
the last 40%, the NACA-001 2 profiles are more favorable. 

Figure 8(a) shows the skin friction coefficient Cf profiles for the upper surface of the CAST7-158 

(solid line) and the NACA-001 2 (dashed line) at y/R = 0.96, 0.90, 0.8, 0.70. The skin friction coefficient is 
a measure of the frictional forces on the surface of the airfoil. In figure 8(a) a smooth variation of Cf as a 

function of x/c up to x/c * 0.4 for the NACA-001 2 and x/c * 0.7 for the CAST7-158 at all four y/R locations 
is shown. The drop in the value of Cf is the indication of shock-induced boundary layer separation and 

possible flow reversal as is shown in the figure for all y/R locations. At y/R = 0.96, 0.90, 0.80, small 
negative values of Cf exist which indicate that the boundary layer is not only separated, but also there is a 

region of reversed flow which gradually weakens as it moves toward the inner region of the blade to 
reattach itself with the flow. 

Figures 8(b) - 8(d) are the displacement thickness 5*, momentum thickness 9* (both normalized 
with respect to chord length), and the shape factor H profiles for the upper surface of the CAST7-158 
(solid line) and the NACA-001 2 (dashed line) at y/R = 0.96, 0.90, 0.80, 0.70. Displacement thickness is a 
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measure of the distance the free stream airflow is displaced outward (i.e., away from the surface) due to 
the decrease of the velocity in the boundary layer. In analogy to displacement thickness, a momentum 
thickness can be defined in accordance with momentum principal which can be established by equating 
the-momentum flow toss due to the wall friction in the boundary layer to the momentum flow in the 
absence of the boundary layer. A value of the shape factor H > 2.5 is a qualitative measure of flow 
separation and reversal. 

In figures 8(b) - 8(d) a smooth variation of 6*. 0*. and H as a function of x/c up to x/c <= 0.4 for the 
NACA-0012 and x/c = 0.7 for the CAST7-158 is shown. An increase in the value of 8*. 0* and H is the 
indication of boundary layer separation and flow reversal and is shown in the figures for all y/R locations. 
Note that in the tip region, the values of 8*. 0* and H at y/R = 0.96 are much larger than values at y/R = 0.70 
due to the higher transonic streamwise velocities. 

Global Flow Features 

Figures 9(a) - 9(d) and 10(a) - 10(d) are the static pressure and Mach contour fields for both airfoils 
at y/R = 0.96, 0.90, 0.80, 0.70. The free stream pressure and Mach contours in these figures are 
normalized to p*, = 1 and M^ = 0.86. These figures show that for both airfoils the shock position is nicely 

captured, and for the NACA-0012 the last 60% of the airfoil's upper surface suffers from flow separation, 
while for the CAST7-158, only the last 30% of the upper surface is separated. The longer separated 
region of the NACA-0012 contributes to larger unsteady pressure oscillations in the aft portion of the flow 
and in the near and far wake, and contributes to stronger noise propagation. 

IV. CONCLUSIONS 

The aerodynamic toads on a onebladed helicopter rotor in forward flight at transonic tip condition 
were calculated. The unsteady, three dimensional, time accurate, compressible Reynolds-averaged thin 
layer N-S equations are solved in an inertial coordinate system on a body conformed curvilinear grid of C-H 
topology. Detailed boundary layer and global numerical comparison of NACA-0012 symmetrical and 
CAST7-158 supercritical airfoils are made under identical forward flight conditions. It was shown that in all 
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spanwise locations the method is capable of resolving the viscous layer, capturing the shock position, and 
predicting the overall flow field behavior. 

The near future plans for the present method is to modify the inefficient C-H grid topology to that 
of either an 0-0 or 0-H topology and to extend the spatial accuracy of the scheme to higher orders to 
resolve the possible attenuation of propagating waves in the far field. 
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y/R 

(CAST7-158) 

Ci 

(NACA-0012) 

Cl 

(CAST7-1 58) 
c d D c d f c d 

(NACA-0012) 
c d D c df c d 

0.96 

1.0602 

0.5552 

0.05871 0.00112 0.05983 

0.03239 0.00104 0.03343 

0.90 

1.0241 

0.5164 

0.05196 0.00111 0.05307 

0.02776 0.00105 0.02881 

0.80 

0.9151 

0.3979 

0.03455 0.00110 0.03565 

0.01625 0.00108 0.01733 

0.70 

0.7910 

0.2605 

0.01852 0.00109 0.01961 

0.00836 0.00111 0.00947 


Table 1 - Lift and drag calculations based on thin layer N-S computations. 
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V„ . WAKE INDUCED VELOCITIES 



Fig.l - Wake modeling by an angle of attack approach 
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Fig. 3 - Partial view of the C-plane of the C-H grid block for CAST7-158 and NACA-0012 airfoils 



Fig. 4 - Grid and boundary conditions for a rotor blade computations 
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PHYSICAL 

DOMAIN 


OUTER COMPUTATIONAL DOMAIN 



Fig. 5 - Physical and transformed computational planes for the C-grid 
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Computed surface pressure distributions for CAST7-158 and NACA-0012 airfoils 
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